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Abstract 

We study the thermodynamics of an exactly solvable model of a self-interacting partially directed 
self-avoiding walk (DSAW) in two dimensions, when a force is applied on one end of the chain. 
The critical force for the unfolding is determined exactly, as a function of the temperature, below 
the G-transition. The transition is second order and characterized by new critical exponents which 
are determined by a careful numerical analysis. The usual polymer critical index v on the critical 
line, and another one, which we call take a non-trivial value which is numerically close to 2/3. 
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I. INTRODUCTION 



The nature of the collapsed phase that a polymer attains in poor solvent conditions is 
still under debate (see e.g. [1, 2] and references therein). Recent experiments on pulling 
of polymers and biopolymers (see e.g. [3-5]) have enhanced the theoretical interest on the 
unfolding transition a collapsed polymers undergoes when subjected to an external force, /, 
applied at its extrema. Until very recently, most of the existing studies on this subject dealt 
with refined version of the mean field studies originally proposed in Ref. [6]. A common 
characteristic of such studies is that, for a self- attracting polymer, they predict a first order 
phase transition in any dimension at a critical force f c (T). At temperatures below the 
0-transition, where the self-attraction prevails, and for an applied force less than f c the 
polymer is in a compact phase. For forces greater than f c the self-attraction is unable to 
maintain the polymer in its compact conformation and the polymer chain is stretched along 
the force direction. However in d — 2 extensive Monte-Carlo simulations [7], performed on a 
self-avoiding walk (SAW) model, suggested that the transition is second order. An exactly 
solvable model, on a lattice of fractal dimension 2, has been analyzed in Ref. [8] and a 
second order transition was found at a critical force f c (T). In Ref. [9] a rationale was given 
for the change in order of the transition as the spatial dimension, d, goes past 2, by means 
of a renormalization group based argument. Within that framework, it was found that, near 
criticality, the projection of the end-to-end distance along the force direction per monomer 
goes like / — f c , near the phase transition, where / is the force and f c is the critical force. 
Numerical uncertainties are too big to critically test this prediction in the self-avoiding walk 
(SAW) model of Refs. [7, 9]. Furthermore it was predicted [8, 9] that above two dimension 
the transition line f c (T) shows a re-entrance at low temperature, i.e. f c (T) increases at 
low T and after reaching a maximum it decreases becoming zero at Tg, the G-transition. 
This behaviour is similar to the one found in theoretical models of pulling of double strand 
DNA where, however, the re-entrance was found at all d [10-12]. Here we will consider a 
simplified polymer model where the chains are represented by partially directed walks, i.e. 
steps with negative projection along the x— axis, (1,0), are forbidden. This model proved 
to be helpful in the past in order to find the phase diagram in the (temperature, fugacity) 
plane for a simplified 9 transition [13-17]. We take advantage of previous contributions and 
generalize the model to the presence of a pulling force along the direction (1,0) (see also 
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Ref. [18]). Surprisingly in this version the critical force as a function of the temperature 
T can be found analytically. With transfer matrix techniques we find that the end-to-end 
distance per monomer goes like (/ — /c) 1 ^ -1 , with £ < 1. With a sophisticated enumeration 
technique [16] we show the correlation critical exponent, u, takes on a non trivial value on 
the critical line, numerically very close to vq = 2/3, the exponent at the G transition. It is 
not clear whether this is an accidental degeneracy or if it can apply also in the undirected 
case, too. For example in the 3d Sierpinski gasket an exact renormalization leads to a non 
trivial /-dependence of v [8]). 

Our work is structured as follows. In Section II, we introduce the model and the basic 
quantities of interest. In Section III, we outline how the transfer matrix can be applied to 
our model, find explicitly the phase diagram (critical line) and give a rough estimate of the 
exponent (. A scaling argument is proposed to suggest that at criticality v — £. In Section 
IV, we review the enumeration technique proposed in Ref. [16], which we use in Section V 
in order to estimate the value of v on the critical line. In Section VI, we critically analyze 
our scaling ansatz and the hypothesis that v — £. Finally, in Section VII we draw our 
conclusions. 

II. THE MODEL 

The model is a DSAW on a two-dimensional square lattice (see Fig. 1), with (non- 
consecutive) nearest-neighbor interactions. A force f , directed along the same axis of the 
walk, is pulling on one end of the DSAW, the other one being fixed. Given a particular 
configuration C the energy is 

E c = -em - fR x , (1) 

where m is the number of interacting pairs and e the energy per pairs, / the modulus of the 
applied force and R x the longitudinal extension of the walk. Then, the canonical partition 
function can be written 

Ql = Ql(/36,/3/)=5>-^, (2) 

c 

where L is the number of the steps of the walk and (3 _1 = T, is the temperature in units of 
the Boltzmann constant. From now on will set e = 1 without loss of generality. From the 
canonical partition function we construct the grancanonical partition function (generating 
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FIG. 1: An example of DSAW configuration. The quantities yi and the force direction are also 
displayed. Thick dashed lines indicate contacts. 

function), 

oo 

g(T,f,z) = Y,Q L z L , (3) 

z being the step fugacity. The (real) singularity closest to the origin, z c (T, /), of the gener- 
ating function Eq. (3) is related to the free energy per monomer as follows 

In Or 

z c (T, f) = lim (4) 

From the singularities of the generating function when / = 0, a complete phase diagram can 
be extracted (see the next Section and also Refs. [14, 15, 17]). In particular, a discontinuity 
is found in the free energy at a value of T = T e , called the 6 temperature. 



III. TRANSFER MATRIX CALCULATIONS AND PHASE DIAGRAM 

Starting from the definition of the generating function, Eq. (3), when / = 0, we observe 
that it can be conveniently rewritten [13, 14] as 

g(TJ,z) = J20 Lx (T,z)exp(PfL x ), (5) 

where Gl x — J2l Ql,l x z l , Ql,l x being the partition function restricted to walks of total 
length L of which L x steps are along the x direction. 

This is useful because now Q Lx can be written in terms of a transfer matrix, T, of 
dimensionality Ly, where L y is the size of our system along the |/-direction[13, 14]. Such a 
transfer matrix T is defined via its actions on the vectors {vi} i=1 £2 , with Vi = (y^, yi+i), Hi 
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being the height of the site in the %— th row which precedes the right-bound horizontal link 
in that column (see Fig. 1), as follows: 



T(vi,v i+1 ) = exp[(3e(mm{\ri\,\r i+1 \})9(-ri ■ r i+1 )] 
x exp[(|r;| + 1) \nz}. 



(6) 



where = — yi and 6(x) is the Heaviside step function. It can be shown that Q(T, f = 
0,z) develops a singularity when A, the largest eigenvalue of T, goes through 1 [19]. This 
means that for large L x 

g Lx cx(\(T,z)) L *. (7) 
Consequently the force-dependent singularity, z c (T, f) occurs when 

A(7»exp(/3/) = l. (8) 

Eq. (8) has a rather deep consequence. In order for the critical fugacity and hence the 
free energy to display a singularity at a non-zero value of the force, i.e. in order for the 
force induced unfolding transition to exist as a thermodynamic transition and not only as 
a crossover, it is necessary that X(T,z), the largest eigenvalue of the transfer matrix when 
there is no force, has itself a discontinuity as z approaches z c (T, f — 0) = z . Otherwise, 
from Eq. (8) it is clear that there can be no such singularity. If there is a transition, then 
we get the following equation for the critical force: 

f c (T) = -T hm lnA(T, z). (9) 

In Eq. (9) the value of A to be put in the right hand side of the equation is the one pertaining 
to the infinite system. A(T, z) for z slightly less than z is plotted in Fig. 3 with a lateral 
size L y up to 40. It is rather clear that a singularity has to be expected at z = zq in the 
infinite size limit. It was indeed shown [14, 15] that for T < T e (T e = 0.8205... in this 
model) there is a singularity of the grand partition function for z = zq = exp(— (3) and for 
this value of the fugacity the biggest eigenvalue is strictly smaller than 1, being: 

A(/3, z = z = exp(-/3)) = \(J3) = j°il±^g) . (i ) 

1 - v^o 
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FIG. 2: Exact critical line as in Eq. (11) together with points corresponding to estimates with the 
transfer matrix calculation, with strip size L y = 20, 40. 

The transition temperature is obtained when X(f3, z = Zq = exp(— /3)) = 1. Consequently, 
the critical line, f c (T), is obtained by putting A(/3) = exp (—/3f c (T)), i.e.: 

l-ex P (-/3/2) 



fc(T) =Tln 



(11) 



_exp(-/3)(l + exp(-/3/2)). 

and is plotted in Fig. 2, where also the results obtained with the transfer matrix with system 
size up to L y = 40 are displayed. In view of Eqs. (8) and (9), we can define a new critical 
exponent ( which characterizes the directed self-avoiding walk. From Eq. (9), if the largest 
eigenvalue approaches its limit value according to the law: 



\(z ) - X(z) ~ (z - z) c 



then one straightforwardly obtains (via Eqs. (8) and (9)) 
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lim ~ (/ - / c (T))VC-i 



(13) 



where (R X (L)) is the average projection of the end-to-end distance of the DSAW along 
the axis (1,0). From Fig. 3 we estimated 1/2 < ( < 1, with ( ~ 0.7 though a precise 
determination is difficult. If £ < 1 in particular the transition is second order. It is widely 
accepted that for d > 2 the transition is first order and so ( = 1. In Ref. [9] a renormalization 
group based argument in d = 2 on the other hand gave for the (undirected) SAW ( — 1/2. 
This argument would apply also to the present case. Given that the transition is second 
order in our model, it is also sensible to look for the value of the critical exponent v (defined 
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FIG. 3: Plot of the largest eigenvalue (T = 0.4 < Tg) versus z. It is apparent that the largest 
eigenvalue approach a limit value as z approaches zq = exp(— f3) (< 1) from below. Thus a 
transition exists in the thermodynamic sense. 

as R g ~ L u , for large number of steps L, where R g is the gyration radius of the L— site 
polymer. In Section VI, using a scaling argument, we shall demonstrate that v = (. 

In the next Section, we shall study the complete canonical partition function, Eq. (2), 
using a powerful method of exact enumeration introduced in [16], that allows us to reach 
large values of L. 

IV. THE METHOD OF ENUMERATION 

As already said, the configurations of the model are directed walks on a two-dimensional 
square lattice with nearest-neighbor interactions. For convenience, we demand that these 
walks end with a horizontal segment. Since the walks are directed in the x-direction we can 
describe these configurations through the distance between two horizontal steps, measured 
in the positive y direction. Thus, we associate to each configuration an A-tuple (r±, r 2 , rjv) 
corresponding to a configuration of total length L = J^li \ r i\ + N. 

The energy due to the nearest-neighbor interactions for each of these configurations is 
(see Eq. (6)) 

JV-l 

U{r u r 2 ,...,r N ) = - ^ min(|rj|, \r i+1 \)9(-ri ■ r i+1 ). (14) 
i=i 

In the following, we assign weights x for steps in the horizontal direction and y for steps 
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in the vertical direction. Then, the canonical partition function is 



Q L (x,y,u) = J2( xe/3f ) 



N 
N=l 



X y L - N uJ U ^-' rN \ 
\n\ + \r 2 \ + ... + \r N \=L-N 

(15) 

where u = exp(/3). 

Now, it is convenient to consider the partition functions Z^ = Z^\x,y,u) for walks of 
total length L + 1 which start with a vertical segment of height r. Then, we have 

L 

Q L+1 (x,y,cv) = ( 16 ) 

r=-L 

(note that Z^ = xQi(x,y,uj)) which satisfies the following recursion relation 

Z { p=xy^h rlL + e^ i (17) 

[ s=-L+\r\+l J 

obtained concatenating these walks. In the Eq. (17), r = —L,...,L, with L = 0,1,2,.... 
Using the symmetry Z^ = Z^ L r \ the Eq. (17) can be written only for non negative r as 

= xfL L + e^ztr-,+ 

I s=0 
L-r-1 } 

+e? f ^ min(r,s) 4-r-i \ ■ ( 18 ) 

s=l J 



Setting x = y = 1 in (18), we obtain the iteration scheme and the free energy J~l{oj) 

PL 



~\nZ^\ The average longitudinal length of the walk (R X (L)} is simply 



(R x (L)) = ^j- ) lnZf. (19) 
Then, we shall proceed as follows: 

1. we calculate the free energy using the iteration scheme proposed in Eq. (18); 

2. using the Eq. (19), we determine how the quantity (R X (L))/L varies against the 
applied force /. 
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FIG. 4: Plot of {R X (L))/L vs. / for various length, L, of the walk. 
V. RESULTS 

The plot of {R X (L))/L vs. / for various values of L is represented on Fig. 4 at T = 0.4 
which is below the O transition occurring at Tg ~ 0.8205.... From Eq. (11) we have 
f c (T = 0.4) ~ 0.764... 

From a careful examination of Fig. 4, we deduce that the quantity (R X (L))/L decreases 
as L v ~ l , where the critical exponent v might depend on the temperature T. In particular 
the data are consistent with v < 1 if / < f c (T) and v — 1 if / > f c {T). In order to find 
more precise values for the critical exponent, we shall proceed along the same lines of Ref. 
[16]. 

An estimation of the critical exponent through the use of the Pade approximants [20] is 
given in Table I. 

Our estimate of the critical exponent, at the critical force, is close to | , the v value at 
the 6 point at / = [15, 17]. 

To get a deeper insight, let us define an L— dependent critical exponent v(L) through the 
formula 

„, n _ ln(^(L)(L + l))-ln(^(L)(L)) 

U[L) ~ ln(L + l)-lnL (20) 

Plotting v{L) versus an estimated correction-to-scaling term a careful extrapolation to L — > 

oo can be performed, determining the critical exponent v for all the values of the force. Let 
us consider three different regimes: 



1. / < fc 



/ 


V 


L 


<fc 


0.501(7) 


< 1900 


= fc 


0.68(8) 


< 1900 


>fc 


1.00000 


< 1300 



TABLE I: Estimates for the critical exponents from a Pade approximants analysis. Note that in 
the f > f c case the error is completely negligible. 

-0.494 p 
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-0.496 - 
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1/L 0.5 

FIG. 5: Plot of v - 1 vs. l/L - 5 , for f < f c and L up to 3000. 

As an example let us consider / = 0.4. We have found that successive estimates for 
the exponent v — 1 with increasing L follow a straight line when plotted against a 
correction-to-scaling term of l/L - 5 (see also the case of Ref. [16] at / = 0). The plot 
is shown in Fig. 5. The extrapolated value for L — > oo gives v — 1 ~ —0.4998, then 
v — \, the exponent typical of a compact phase. 

2. / > U 

As before, we have plotted the exponent v — 1 against a well-defined correction-to- 
scaling term. Now, this term is order of l/L. Fig. 6 shows the / = 1.0 case as typical 
example. Now, the extrapolated value gives v — 1 ~ 3.0 x 10~ 6 , then v = 1 within the 
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FIG. 6: Plot of v - 1 vs. 1/L, for f > f c and L up to 3000. 
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FIG. 7: Plot of v - 1 vs. 1/L - 28 , for / = / c and L up to 3000. 
numerical precision. 
3. / = U 

Now, the correction-to-scaling term is of 1/L 28 (see Fig. 7) and v — 1 ~ —0.3336, 
which implies z/ ~ |. 

Thus we find that the value of v for / = / C (T) is equal (in the limit of numerical precision) 
to the value |, which corresponds to that for / = at T = Tg [15, 17]. This is a non trivial 
result. In particular, we have to expect that along all the critical line / = f c (T), v takes 
the value |. Moreover, as pointed out in [17], the correction-to-scaling term when / = 
and T = T e is of 1/L 1 / 3 . In our case, we found that this correction increases to 1/L - 28 , 
for T = 0.4. Within the numerical errors, this implies that the correction to scaling scaling 
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<fc 


1/2 


= fc 


2/3 


>fc 


1 



TABLE II: Estimates for the critical exponents from an extrapolation to L — > oo, obtained plotting 
the L-dependent critical exponent v(L), Eq. (20), versus an estimated correction-to-scaling term. 

term depends on force as well as on temperature and that the size of our system has to 
increase in order to find the right critical exponent. 

Finally, in Table II we have summarized the above results. 

In the next Section, we shall introduce a scaling theory that rationalizes for what we have 
found on the critical behaviour of the average horizontal end-to-end displacement (R X (L)). 



VI. SCALING THEORY 

Our previous results suggest the following scaling ansatz (see also Ref. [7]): 

(R X (L)) = L v Q(AfL+) (21) 

where Af = f — f c (T). The scaling function <J>(x) must have the following behaviour: 

' x (i-v)/i> if a; -> +oo 



~ <( $ if x -> , 

\ x i(i/2-v)W if x _^ _ QO 



(22) 
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with $o a not-zero constant value. Then, the quantity (R X (L)) obeys to the equations: 



LAf^ 1 -^ if Af > 0, if > 



(R X (L)) ~ { 



if A/ = 



k LV2| A/ |(i/2-.)/v if Af < 0, ^ < 

in agreement with the results found in the previous Section. 

Now, let us observe that the free energy contribution to the singular part is 

AF = (R x (L))f - (R x (L))f c = (R x (L))Af 
= AfL v <5>(AfL^) 



(23) 



(24) 



where we have used Eq. (21), / is the applied force and $(x) = x&(x). Since AF is a 
contribution to the total free energy (not a free energy density), we expect it depends only 
on the "dimensionless" combination of the scaling fields Af and L with the appropriate 
exponents. This implies that v — xjj. 

Comparing Eqs. (13) and (23), we deduce that ( = v = 2/3. Then, Eq. (21) becomes 
(R X (L)) = L 2 / 3 $(A/L 2 / 3 ). 

To test this prediction, we have plotted in Fig. 8 the function L' 2 ^ (R X (L)) versus / — / c , 
where f c is again determined from the exact formula, Eq. (11). It is evident that, apart 
from obvious finite size scaling corrections, our ansatz is justified. Then, we have derived 
the scaling function $(x). The final result is shown in Fig. 9. 



VII. CONCLUSIONS 



In this work we have analyzed a model of self-avoiding partially directed chains on a 
square lattice, with a force pulling along one of the two lattice directions. The model is 
simple enough to yield the exact form of the temperature dependent critical force, / C (T) (see 
Eq. (11)). However, the critical indices of the unfolding transition, which is second order, 
are not trivial. The transition is characterized by two exponents, the usual correlation 
length critical exponent u, and one which we called (. In particular, the exponent v at 
f — fc ( see Eq. (23)) is different both from 1/2, the collapsed polymer value, and from 1, 
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FIG. 8: Plot of (R X (L)) scaled by L 2 / 3 versus / — f c . Finite size scaling corrections to the critical 
force are evident. 
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FIG. 9: Successive estimates for the scaling function are shown for increasing lengths of the 
walk. 

the extended polymer value (see again Eq. (23)). The ( exponent characterize the singular 
behaviour of the chain elongation per monomer along the force direction as the critical 
force is approached from above (see Eq. (13)). Through a powerful enumeration technique 
taken from the literature [16], coupled with a finite size scaling to extrapolate our results to 
infinitely long chains, we find that v is very close to 2/3. A scaling analysis also suggests 
that v — £ at least within our numerical precision. Further investigations are required to 
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extend our results to the undirected SAW case. 
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